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We discuss the soft contribution to the elastic pion form 
factor with the mass evolution from current to constituent 
quark being taken into account using a light-front Bethe- 
Salpeter (LFBS) model, which is a light-front quark model 
(LFQM) with a running mass. It is shown that partial con- 
servation of the axial- vector current (PCAC) is satisfied with 
a running quark mass. We examine the sensitivity of the pion 
form factor using two different functional forms of the quark 
propagator. The Ball-Chiu ansatz is used to maintain local 
gauge invariance of the quark-photon vertex. The extension 
of our model to the hard contribution is also discussed. 

PACS numbers: 12.39.Ki, 13.40.Gp, 14.40.Aq, 14.40.Lb 



I. INTRODUCTION 

The pion electromagnetic (EM) form factor is of great 
interest for the study of Quantum Chromodynamics 
(QCD). At low momentum transfers (Q 2 ) nonperturba- 
tive QCD (NPQCD) dominates, while at large Q 2 per- 
turbative QCD (PQCD) can be used to calculate the 
asymptotic form factor; and the transition from NPQCD 
to PQCD has long been of interest 0. Since PQCD 
can only be used for Q 2 > 1 GeV 2 , light-front (LF) 
quantization methods may be most useful [Q. In the 
early work on the pion form factor || using the light- 
front Bethe-Salpeter equation (LFBSE) {§ the NPQCD 
(soft) part was separated from the PQCD (hard) part; 
and it was shown in a model calculation that the tran- 
sition from NPQCD to PQCD is expected in the region 
5.0 < Q 2 < 15.0 GeV 2 . This work was extended || to 
include the Sudakov form factor S, anomalous quark 
magnetic dipole moments, and a simple model for the 
running quark mass, m(Q 2 ). Improved calculations of 
the pion EM form factor are motivated by the recent new 
and upgraded data (up to Q 2 =1.6 GeV) from the Jeffer- 
son Laboratory(JLab) Q]. While this Q 2 range may be 
too low to determine the transition to the PQCD region, 
these data are useful for studying NPQCD theoretical 
approaches. In the present work we restrict ourselves to 
the soft NPQCD part with models for the running mass 
that ensure gauge invariance and consistency with the 
partially conserved axial current (PCAC) relation H 
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quark condensate. 

Since we are only considering the soft part of the pion 
form factor here, the LFBS amplitude can be modeled 
by a light-front wave function || based on LF Hamilto- 
nian dynamics, such as the light- front constituent quark 
model(LFCQM) ]lO|-[T^|, but an essential ingredient is 
the use of a running quark mass, which is the main sub- 
ject of the present paper. In LF quantization, a pos- 
sible connection is anticipated between the constituent 
quark model(CQM) and QCD due to the rational energy- 
momentum dispersion relation that leads to a relatively 
simple vacuum structure. There is no spontaneous cre- 
ation of massive fermions in the LF quantized vacuum. 
Thus, one can immediately obtain a constituent-type pic- 
ture, in which all partons in a hadronic state are con- 
nected directly to the hadron instead of being simply 
disconnected excitations (or vacuum fluctuations) in a 
complicated medium. In particular, a systematic pro- 
gram has been laid out in the calculation of the space- 
like EM form factor of pseudoscalar mesons because only 
parton-number-conserving Fock state (valence) contribu- 
tion is needed when the "good" components of the cur- 
rent, J + (= J° + J 3 ) and J± = [Jx, J y ), are used in the 
Drell-Yan-West(<7 + =0) frame fll3|Ji- The new data from 
the JLab |Q seem to be in good agreement with the pre- 
vious CQM result based on a QCD-motivated linear 
confining potential. However, a possible realization of 
chiral symmetry breaking in the LF vacuum is an under- 
developed aspect of LF quantization. 

In contrast to quark models or LFCQM which use a 
phenomenological constant constituent quark mass, an 
approach based on QCD quantum field theory is the 
Bethe-Salpeter (BS) equation in conjuncton with Dyson- 
Schwinger (DS) equations for the quark propagators, 
gluon propagator and vertices. We note an important 
result of recent DS calculations, in which the effective 
running mass, m(p 2 ), is calculated juj. In these DS cal- 
culations the parameters for model gluon propagators are 
fixed by fitting the quark condensate, the mixed conden- 
sate and even the form of the nonlocal condensate [ p^| . 
It is particularly interesting to note that the DS quark 
propagator with running masses that decrease with re- 
spect to p 2 faster than quark models gives properties of 



where muA=mQ (current mass) in spacelike p 2 =v 2 
region(f is the renormlization point) and (qq)t v ) is the 



the rho (l6| and pion ||17[ that are in agreement with 
experiment. In all of these DS calculations it has been 
found that the effective quark mass drops very rapidly 
with increasing Q 2 . Since the covariant BS/DS approach 
and the LF approach are not same but complementary, 
it may be necessary to examine if the previous LFCQM 
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result |12] is intact even if the quark mass evolves as 
rapid as BS/DS approach found. This is a strong mo- 
tivation for reformulating the LFCQM with a running 
quark mass, which we do in the present work. 

In the present work, we analyze the effect of the mass 
evolution (from constituent to current quark mass) on 
the elastic pion form factor at low and intermediate Q 2 . 
A correct theoretical approach to such a study is the LF- 
BSE coupled to LFDS equatons. From the LFDS equa- 
tions the running quark mass is obtained from the dressed 
quark propagator. Although the LFDS equation has re- 
cently been developed [lSy , only simple model solutions 
are available, and here we use a parameterization of the 
running quark mass that is consistent with known ob- 
servations. While the asymptotic behavior of the run- 
ning mass might require the crossing symmetry (under 
Q 2 «-> — Q 2 ) at high momentum Q 2 analogous to that of 
pion form factor, there is no clue yet for the small mo- 
mentum behavior in timclikc region. Thus, we present 
the two different forms of mass evolution function; one is 
crossing asymmetric and the other is crossing symmetric. 
We then compare the results for the two cases. It should 
be noted that the recent LFDS results [|l8| also show a 
rapid decrease of effective quark mass with momentum 
that puts the use of quark models for calculating any but 
static properties in question. In the present work we use 
forms for the quark mass evolution that are consistent 
with conventional quark models for calculating hadronic 
properties at momentum transfers less than about 800 
MeV. 

The paper is organized as follows: In Sec. II, we review 
the formulation that underlies a description of the elastic 
pion form factor within a modeling of QCD through the 
LFBSE, which is a LFQM with a running mass. In Sec. 
Ill, we formulate the running quark mass in LF frame- 
work and discuss the local gauge invariance at the quark- 
photon vertex, i.e. Ward-Takahashi identity fisf] , due to 
the momentum dependent quark propagator. In Sec. IV, 
we analyze the running mass effect on the pion form fac- 
tor, charge radius, and decay constant obtained from the 
previous CQM calculation. We also show that our 
model with the running quark mass is consistent with 
the PC AC relation given by Eq. (0). This simply means 
that we obtain a quark condensate consistent with the 
phenomenological value given in Eq. ([l]) , which is true in 
the DS calculations |T4| , |l5| ]. A conclusion follows in Sec. 
V. 



II. MODEL DESCRIPTION 

The LFBSE introduced in Refs. §,§ has the form 

fcc(x,kr-,y, Ll) 



*(x,k ± ) = J[dy}[d 2 l ± ] 



where IC g is the gluon exchange kernel and /C c , the differ- 
ence between the complete BS kernel and /C g , includes all 
confining effects. These kernels were obtained from the 
relativistic string and LF perturbative QCD (PQCD), re- 
spectively. Since NPQCD has not yet provided a form 
for the confining kernel, in Ref. Q the problem of solving 
for the complete BS amplitude, ^(x, kj_), was avoided by 
using a model for the soft amplitude. I.e., the soft BS 
amplitude, ^(a;, kj_), can be considered to be the solu- 
tion of the equation 

* s (x,k x )= [[dy][d 2 l ± ]lC c (x,k ± ;y,l ± )* s (y,l ± ). (3) 

Iterating Eq. (^) by inserting ty s for ty, one obtains the 
approximate form 

*(a:,kx) «**(*, k x ) 

+ f[dy][d 2 l x }]C s (x,k ± ;y,lxW(y,lx). (4) 

This BS amplitude contains both soft and hard ingre- 
dients needed to take care of momentum transfer for 
all Q 2 therefore is correctly characterized as including 
both confinement and asymptotic features of a compos- 
ite quark system. This approach to the pion form factor 
has been shown j|] to be in good agreement with the 
direct BS calculation and to converge rapidly. The 
extension of Eq. (|^) to the non-wave- function vertex in 
the particle-number-nonconserving Fock state contribu- 
tion has recently been discussed in Ref. ^3] . 

One nice feature of this approach is that one deter- 
mines the soft part and the hard part separately, so that 
one can determine the transition from soft to hard QCD 
within this LFBS approach. This will be the subject 
of our future work. For the present work of low- and 
medium-Q 2 , we only consider the confining part of the 
BS amplitude. We thus may be able to use the LFCQM, 
which has inlcuded many important properties of the Q 2 
range that we focus in this work. Therefore, in the rest 
of this work we use the LF wave function for the LFBS 
amplitude. 
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FIG. 1. The valence diagram with dressed quark propaga- 
tors for the pion EM form factor calculation in q + = frame 
where p q = P — k, p' q — p q — q — P' — k and p q = —k. 

The elastic pion form factor is related to pion EM cur- 
rent by the following equation 



*(l/,U), (2) 



W(0)|P) - {P' + PYF 7 ,{Q 2 ). 



(5) 
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As usual, our calculation will be carried out using the 
q + = frame where q 2 = (P-P') 2 = q + q~ -q 2 ± = -Q 2 , 
i.e. Q 2 > is spacelike momentum transfer. 

The matrix element of the current given by Eq. (|^) 
can be expressed as a convolution integral in terms of LF 
wave function, ty s (x, kj_) as shown in Fig. [IJ 
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dx / d 2 k ± n*x q -(*y±) 
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where Pg=p' q + = {\ — x)P + and k'x=kx — a;q± in the 
initial pion rest frame, Px=0. The helicity of the 
quark(antiquark) is denoted as X q igy Since the matrix 
element of the current in Eq. (^|) is symmetric under the 
exchange of q and q, we do not explicitly write the con- 
tribution of photon-antiquark interaction diagram as well 
as the charge factor. 

Our LF wave function ^> s in Eq. (^J) is given by 
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n q X, (X, kj.) = y -^*(X, k X )^A,A 4 (X, k X ), (7) 

where $ and TZ are the radial and relativistic spin-orbit 
wave functions, respectively. Our radial wave function 
is given by the gaussian trial function for the variational 
principle to the QCD-motivated effective LF Hamiltonian 



$(k 2 ) = ( 



1 



) 1 / 2 exp(-k 2 /2/3 2 ), 



(8) 



where k=(fc z ,kx) is three vector and <&(k 2 ) satisfies 
fd 3 k\$(k 2 )\ 2 = 1. The LF variable (x,k ± ) is introduced 
in Eq. (ra) by the definition of the longitudinal momentum 
k z via k z =(x - 1/2)M with M§=(m 2 + k 2 ± )/x(l - x). 
If the quark mass depends on x and kxj the Jacobian of 
the variable transformation (k z , kx)— >(x, kx) in Eq. (Q) 
is obtained as 



dx 



Mi 



(2a; — l)m(a;, kx) dm(x, kx) 



(9) 



4x(l - x) 2x(l - x)M <9x 
The spin-orbit wave function for a pseudoscalar meson 



(J PC = + ) is obtained 10 12] by the interaction inde- 



pendent Melosh transformation as follows: 

u{jpq,\q)lbV{Pq, Xg) 
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III. QUARK MASS EVOLUTION AND LOCAL 
GAUGE INVARIANCE 

The solution of the DSE for the renormalized dressed- 
quark propagator takes the form in Minkowski space 



s( P y 



A(j?)i>-B{p 2 



(11) 



where the quark mass evolution function m{p 2 ) is defined 
as m(p 2 ) = B(p 2 )/A(p 2 ). Also, the gauge invariancc 
requires that the quark-photon vertex T^ 1 given by Eq. (|^) 
satisfy the vector Ward-Takahashi identity(WTI) |g (i.e. 
current conservation) 



s( P y 



(12) 



where q = p — p' . At zero momentum transfer q=0, 
the quark-photon vertex is also specified by the differen- 
tial Ward identity (i.e. charge conservation) T tl (p;Q) = 
8S(p)^ 1 1 'dp p. The bare quark-photon vertex, T M = 7 M , 



which is usually used in LFCQM |10|-|12[, is inadequate 
when the quark propagator has momentum-dependent 
dressing because it violates WTI. As used in many DSE 
studies of EM interactions [^l|jl4| , we take the Ball- 
Chiu(BC) ansatz |22] for the quark-photon vertex 
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,^ A(p' 2 )~A( P 2 ) 



p/2 _ p2 



+ A(p' 2 ) + A(p 2 ) r _ {p + pr B(p< 2 



B(p 2 



p/2 _ p2 



(13) 



Here, we introduce two algebraic parametrizations of the 
quark running mass, i.e. crossing asymmetric (CA) and 
crossing symmetric(CS) mass functions proportional to 
p 2 and p , respectively. For the CA mass evolution func- 
tion, we take 



m(p 2 ) = m + (m c - m ) 



l + cxp(- M 2 /A 2 ) 
l + cxp[(-p 2 -fj, 2 )/X 2 Y 



(14) 



where Too and m c are the current and constituent quark 
masses, respectively. The parameters fj, and A are used 
to adjust the shape of the mass evolution. Similarly, the 
following form of the mass evolution function is used for 
the CS case 



m(p 4 ) = to + (to c - to ) 



1 + cxp(- M 4 /A 4 ) 
l + exp[(p4_ M 4)/ A 4]' 



(15) 



where we simply replace —p 2 and /i 2 (A 2 ) in Eq. ( |l4| ) with 
p 4 and /i 4 (A 4 ), respectively. In our case, we set A(p 2 )=l, 
B(p 2 )=m(p 2 ) for CA and B{p 2 )=m{p i ) for the CS case, 
respectively. While Eqs. (|l4|) and (|l^) are phenomeno- 
logical, the results of our running mass in spacelike mo- 
mentum region (— p 2 > 0) yield a generic picture of the 
quark mass evolution from the low energy limit of the 
constituent quark mass to the high energy limit of the 
current quark mass. For comparison, we use in Fig. ^ two 
different parameter sets for each mass evolution function, 
i.e., (m 2 ,A 2 ) = (0.9,0.2) [Set 1] and (0.5,0.2) [Set 2] (in 
unit of GeV 2 ) for m{p 2 ) and (/i 2 , A 2 ) = (0.95,0.63) [Set 
1] and (0.28, 0.55) [Set 2] (in unit of GeV 2 ) for m(p 4 ), re- 
spectively (See Fig. || for the line code). The current 
and constituent quark masses used are toq = 5 MeV 
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and m c = 220 MeV, respectively. Simulating the con- 
stituent picture at small momentum region, we have cho- 
sen these particular sets of parameters, [Set 1] and [Set 2] 
for each mass function, to keep the constituent mass up 
to (— p 2 ) ~ 1 and 0.5 GeV 2 , respectively, before it drops 
exponentially. 
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FIG. 2. Quark mass evolution in spacelike momentum re- 
gion, — p 2 > 0. 

In order to express the four momentum p 2 in terms of 
LF variables (a;, k^), we use the on-mass shell condition, 
p 2 = m 2 {p 2 ). It implies zero binding energy of a mock 
meson, i.e. P~ = Pq + Pq where -P~(= P° — P 3 ) and 
PqiPq) are the LF energies of the mock meson and the 
quark(antiquark), respectively. It leads to the following 
identity for the pion case (m q — rriq) 

p 2 = x{l-x)M 2 -ki. (16) 

For the mock meson mass M, we take the average 
value (so called spin-averaged meson mass) of 7r(m w ) 
and p(jn p ) masses with appropriate weighting factors 
from the spin degrees of freedom, i.e. M—(m n + 
3TOp) eX p/4=612 MeV, which is consistent with quark 
model calculations with typical constituent quark masses 
and been used in some constituent LFQM calcula- 
tions jllj]. Note that the spacelike form factor F- K {Q 2 ) 
is weakly dependent on M, i.e. the numerical result with 
the mock meson mass M = 612 MeV is not much dif- 
ferent from that with the physical pion mass = 140 
MeV, which would be used in the LFBS-LFDS approach. 
Using Eq. (|l|), we can now express the mass evolution 
functions m(p 2 ) and m(p ) in terms of LF variables x 
and k^. 

Using the good component of the current, J + , the pion 
EM form factor in Eq. (|5|) is obtained as 



Ftt{Q 2 ) =N * J ^ 2 kxy^y^$/(^k'x)^(a:,k ± ) 

/ f k± ■ k'± + m k mw m k Am k (2M 2 + q^J \ 
X \ x(l-x)M M^ M M^Ak 2 ± y 

(17) 

where is the normalization constant and Amj = 
m(x,k'^) — ni{x,k^)~ra k i — m k and Ak^ = k'^ — 
k^. The other primed terms, k' z and Mq, are ob- 
tained from k z (x,k± — > k'^) and Mq(x,\l± — > k'^), 
respectively. The terms in the curly bracket are ob- 
tained from the trace of spin-orbit wave function, i.e. 

J2 n H^/ \/l^) r Bc( u / \fp^) n - Note that thc normal- 
ization constant at Q 2 — is exactly one in chiral 
limit (m = M = 0). 

We also obtain the quark condensate from the PCAC 
relation given by Eq. ([I]) as follows: 

,_ > /ttV^Vc f dxd 2 k± ldk z / 

(18) 

where N c =3 is the color factor. This is equivalent to 
the expression in DS models [jlj] for A(p 2 ) = 1. Thc 
quark condensate is normally evaluated at the spacelike 
momentum scale p ~ 1 GeV(corresponding to the renor- 
malization point v |14|| ) where m k — > too (see Fig. ||). 

IV. NUMERICAL RESULTS 

In our numerical calculations, we use the model 
parameters (to c , /3)=(0. 22, 0.3659) [GeV] obtained in 
Ref. (l^] for the linear confining potential model where 
the charge radius(r 2 =— &(1Ft T {Q 2 ) / dQ 2 \Q2 =Q ) and decay 
constant((0|g7 AI 75(7|P) = ifnP 11 ) of the pion were pre- 
dicted as r 2 =0.425 [fm 2 ](Exp. = 0.432±0.016 ||) and 
7^=130 MeV (Exp. = 131 MeV §|), respectively. The 
change of the charge radius and decay constant from 
the CQM result due to the running mass formulae are 
within 2% and 5%, respectively. For the calculation of 
the quark condensate, we obtain, for example, 7^= 127 
MeV from the [Set 1] of the CA mass function. Con- 
sequently, we obtain from Eqs. (P and ([l8]) the pion 
mass and the quark condensate as m^—WA MeV and 
~(qq) — (0.3 GeV) 3 while the experimental values of m w 
and ~{qq) are 140 MeV and (0.236 GeV) 3 gf|, respec- 
tively. This shows the PCAC relation is reasonably well 
satisfied in LFQM with our mass evolution function. The 
results from other parameter sets are not much different 
from the above [Set 1] of CA case. 

In Fig. ||, we show our results of the pion EM form 
factor for small Q 2 region using the [Set 2] with BC ver- 
tex for both CA and CS mass evolution functions and 
compare with the experimental data p3]| as well as the 
CQM result given by Ref. |12] |. The line code is in the 
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figure. The small momentum(Q 2 ) behavior of the form 
factor with CA and CS mass functions are not only close 
to each other but also close to the experimental data |||] 
as well as the CQM result fl§]. The results of [Set 1] for 
both CA and CS mass functions are even closer to the 
CQM result. 



used for the dressed quark-photon vertex. We showed 
that the PCAC relation in Eq. (0) is satisfied in LFQM 
with our running mass formulae. The CQM result is not 
affected too much by the quark mass evolution for the 
small momentum transfer rcg ion up to Q 2 =l~2 GeV 2 . 
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□ Amendolia et al. (1986) 
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Bebek et al. (1978) 
• Jefferson Lab (2000) 
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FIG. 3. Pion EM form factor for small Q 2 region using the 
[Set 2] for both CA and CS mass formulae compared with the 
experimental data |?3j and the CQM result p2[ . 

We also show in Fig. [| our results of the form factor 
for the intermediate Q region for CA [Fig. |(a)] and 
CS [Fig. [|(b)] mass functions compared with the exper- 
imental data §J26| as well as the CQM result @. The 
line code is given in each figure. As one can see from 
Fig. H (1) the difference between the bare vertex and 
BC ansatz indicates the breakdown of the local gauge in- 
variance from the usage of the bare vertex, (2) the [Set 
2] for both CA and CS mass functions show larger devi- 
ation from the CQM result than the [Set 1] case for the 
momentum transfer Q 2 ^2 GeV 2 and above region, (3) 
the results with BC vertex fall off faster (at around Q 2 —2 
GeV 2 ) than the CQM result does, (4) the mass evolution 
effects from [Set 1] for both CA and CS cases are not 
much different from the constituent result up to Q 2 —8 
GeV 2 , and (5) the CA mass evolution function is more 
sensitive to the variation of the momentum dependence 
than the CS mass evoluton function. 



V. CONCLUSION 

In conclusion, we have reexamined the soft contribu- 
tion to the pion elastic form factor in the framework of 
the LFBS with a Q 2 -dependent quark mass that could be 
obtained from a LFDSE. This is equivalent to the LFQM 
with a running quark mass. The Ball-Chiu ansatz was 
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" Bebek et al. (1978) 

* Jefferson Lab (2000) 
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FIG. 4. Pion EM form factor:(a) Crossing asymmetry(CA) 
and (b) Crossing symmetry(CS) mass functions compared 
with the experimental data |7jj2d] as well as the CQM re- 
sult @. 

However, the form factor is sensitive to the shape of 
the mass evolution for the intermediate ranges, e.g. the 
[Set 1] for both CA and CS cases are not much different 
from the LFQM result but the [Set 2] show a sizable ef- 
fect on the soft pion elastic form factor for Q 2 ~ 2 GeV 2 
and above. It may be interesting to observe from [Set 
2] that the form factor may distinguish the mass evolu- 
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tion respecting CS from the one not respecting CS even 
though they show similar momentum dependent behavior 
as shown in Fig. 0. While our calculation employed a phe- 
nomenological mass evolution functions and the quark- 
photon vertex modification, the qualitative feature pre- 
sented in this work might not be significantly modified 
even if one were to use the more realistic solutions ob- 
tained from the first principle. However, it may be inter- 
esting to check further whether the soft par t would fall 
more rapidly if the present results of DS (l4| and LFDS 
p8| models were used. The "hard-scattering" contribu- 
tion to the pion form factor, i.e. the second term in 
Eq. (Q), is under investigation. 
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